INFINITE PROPAGATION SPEED FOR WAVE SOLUTIONS ON 
SOME P.C.F. FRACTALS 



YIN-TAT LEE 



Abstract. The finite difli'erence method for the wave equation on p.c.f. frac- 
tals suggests that the propagation speed of the wave equation may be infinite. 
We prove this is indeed true if the heat kernel satisfies a sub-Gaussian lower 
bound. Furthermore, we provide a sub-Gaussian upper bound for the solution 
of the wave equation given the heat kernel sub-Gaussian upper bound. 



1. Introduction 

In [3], Dalrymple, Strichartz, and Vinson pointed out that there is no maximum 
propagation speed on the Sierpinski Gasket (SG) because of a scahng property of 
SG. In other words, there is no C such that for all x and t > 0, the fundamental 
solution of the wave equation at point x and time t is supported in Bct{x). However, 
it does not rule out the possibility that the fundamental solution is supported in 
Bf(^i){x) for some continuous function / such that /(O) = 0. 

In this paper, we first provide an error analysis for the finite difference method 
on p.c.f. fractals with regular harmonic structure. Let u be a solution of the wave 
equation on the fractal K, and let Um be the solution on the level m approximation 
Vm of K. In Theorem ini we show that Umix^t) ~ u(x,hmt) where hm is a time 
renormalization factor. Interestingly, the hm decreases faster than the grid size 
does as m increases for most of p.c.f. fractals. It means that the propagation speed 
of Um increases as m increases. Although the result will not be used in the later 
proof, it gives the heuristic reason why the infinite speed holds. 

In Theorem [51 we prove the infinite propagation speed. If the initial position 
is zero and initial velocity is positive, then u attains positive values for all points 
X € K within arbitrary small time period. The proof uses a heat kernel lower bound 
and a relation of heat equation and wave equation. In Theorem ll2l we prove a off- 
diagonal upper bound for the solution of wave equation using a complex time heat 
kernel upper bound. This upper bound is also sub-Gaussian. 



2. Preliminaries 

At first, we define briefly the notations and concepts introduced by Jun Kigami[TU]. 
An iteration function system (IFS) is a finite set of contraction mappings {Fi}fLi 
on a complete metric space. An IFS fractal K is the unique compact set such that 
K — [j^^^FkK. A connected IFS fractal is called post critical finite (p.c.f.) if 
there is a finite set Vq such that FjK H FkK C FjVq (1 FkVo ior j ^ k. For a word 
Lu ~ 1112 . . .im, we define Fi^ — Fi-^ o Fi^ o ■ ■ ■ o Fi^ . For example, the interval [0, 1] 
is the unique IFS fractal generated by mappings {Fi{x) — ^x,F2{x) = ^x + ^} 
and the corresponding Vq is {0, 1}. 

1 
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Now we define a sequence of increasing finite graplis Ti to approximate K. Let 
To be tfie complete grapli of the finite set Vq. For i > 0, we define Fi = {Vi, Ei) 
where Vi = yj^^-^FkVi-i and 

Ei = {{FjX,Fjy) eViX Vi\ {x,y) e Ei_i}. 

Wc define V, = U^j^Vfe. Using [0, 1] as an example, the corresponding Vn is € 
{0, • • • , 2"}}, r„ is the simple path with 2" + 1 vertices and K is the set of dyadic 
numbers {^|n e N and < o < 2"}. 

For any finite set V, a non-negative symmetric bi-lincar form 5 on is called 
a Dirichlet form if £{u,u) = for all constant functions m on y and £{u,u) > 
£{[u], [u]) for any function uonV where [u] = min{max{u, 0}, 1}. For V C V, we 
can induce a Dirichlet form Sy' from by 

£v'{u,u) = inf £v{v,v). 

v\yi=u 

For any IFS fractal K, a sequence of Dirichlet form {£m} on {Vm} is called com- 
patible if Em is induced from £m+i for all m. If this sequence satisfies the equation 

AT 



£m+i{u,v) = ^£rn{uo Fi,V O Fi) 



i=l 

for some number r, > 0, we call it a self-similar sequence and it is said to be regular 
if n < 1. 

If £m are compatible, £m{u\vrr^ ^ Mv^) is increasing. For any function u on K, 
we define energy £{u,u) as limm^oo £m{u\v^,u\vrn) and dom£ = {u G C{K) : 
£{u,u) < oo}. It is known that dom£/constants is a Hilbert space. For any 
function u on Vm , the harmonic extension of u is the unique continuous function u 
on K minimizing the energy £{u, u). We define ip'^ix) to be the harmonic extension 
of the delta function S^p on Vm- 



For example, we can define a regular self-similar sequence on [0, 1] by 

2"'-l 

£m{u,u) = 



fe=0 



1 

2m 



1 

2m- 



The corresponding energy £{u,u) = \u'{x)\'^dx for u £ C^{[0,1]). The corre- 
sponding harmonic extension is linear interpolation on Vm', ip^ix) is a triangular 
function at point p with width 2~™+^ and dom£ = H^. 
We define the resistance metric on K by 

R{x,y) = msiK{£{v,v)~^ : v{x) = 1, v{y) = 0}. 

It is known that K is compact under resistance metric, in particular, 

(2.1) \u{x) -u{y)f <C£{u,u). 

where C = sup^ y^j^ R{x, y) < oo. 

Next, we define a self-similar probability measure on ii" by 



N 
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for some /i; G (0,1) such that ^ /ii = 1. For u G doni£, the Laplacian of u 
corresponding to the self-similar /i is defined by the weak formulation: A^w = / if 
f £ and 



£{u,v) = — J fv dfi 

for all V G domS that vanish on the boundary Vq. If A^u is continuous, we have a 
pointwise formula for A^u: 

Auu{x) = lini fj.:^^^Hmu{x) 
where /im,p = / i^p^ix) dx and H„i is the self-adjoint matrix such that 

£m{u,v) ^ -'{u,H„iV). 



3. Existence of solutions 

Let i? be a finite subset of Vt. For u G C{K) vanishing on B, we define the 
Laplacian A^^bu = / if ^J,^^H„iu{x) converges uniformly to a continuous function 
/ on K\B. A.2] 

The wave equation with boundary set B, initial position / and initial velocity g 
is defined by 



{x G K\B, t G R) 
{x^K) 
(xeK) 
{x£ B, i G M) 

where the time derivative uu is in the classical sense. For convenience, we write A 
instead of A^_b- The condition B = $ corresponds to Neumann boundary condition 
and B = Vq corresponds Dirichlet boundary condition. In this paper, we use C as 
a generic constant which depends only on the fractal. Since most of the proofs in 
this and next sections need extra care for the case _B = 0, we omit the proofs for 
that case in these two sections. 

By [121 A.2], we have a set of orthogonal eigenvectors of —A with 

corresponding increasing eigenvalues {A„}„>]^ such that Hvnlb = 1 and {<Pri}„>i 
spans domA. By the assumption i? 7^ 0, we have A > 0. 

Lemma 1. dom£ — {J2'^n'Pn ■ J2'^n'^n < 00} and ||m||oo < C£(w)^/^. 

Proof. Let u — ^a„iy9„ such that M = ^a^A„ < 00. Let Um be the partial sum 
of u. Using £{um) — SnLi '^nAn < M and p.l|) , we have 

\uM~u„^iy)\ < C M^/^. 

Combining with J it^„ < M/Ai, we get 

lk™||oo < (C + A7'/') M^/\ 

Hence, Um is uniform bounded. Again, by (j2.ip . Um is equicontinuous. So, by the 
Arzela-Ascoli theorem, u G C{K). Since dom£/constants is a Hilbert space, we 
get £{u) = limf (wm) = M < 00. Hence u G dom£. 

The converse follows from £{Y^ anfn)=J2(^n^n- O 



{utt{x,t) = Af^^Bu{x,t) 

u{x,0) = fix) 

ut{x,0) = g{x) 

u(x,t) = 
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For domA, we do not have a similar description using the original definition. 
So, we extend the domain of A to {'^an^n ■ ^o-n^n < ^} by the identity 
A(^a„(/3„) = —'^(inXn^n, which converges in L^. Since ^a„(/9„ converges in 
L°° and ipnls — 0, the boundary condition is satisfied for u £ doniA. 

Let the initial position be / = ^ a„(y9„ and the initial velocity be g = ^ 
We define the formal solution by 



, . ^ — T , rr — , ^ — T sin(\/ A^i) 

(3.2) U = 2_^anCOS[y'Xnt)(pn+ 2_^Pn 7= 

" 'An 



It is standard to prove the formal solution is a weak solution under some condition 
on / and g. In [9], Hu discussed wave solutions for the Frechet derivatives. However, 
in order to complete the error analysis using the finite difference method, we need 
to prove that the classical solution exists. 

Theorem 2. /// G rfomA with Af G dom£ and g G domA, then the solution u of 
the wave equation exists. 



Proof. Let u be the weak solution defined by (|3.2p . Formally, we have 
utt = -^a„A„cos(\/A^t)(^„ + /3n\/A^sin(\/A^t)y)n. 

n n 

Fix xo £ K and 7„ = Q;nA„sgn(iy9„(a;o))- Since Af G domf , we have 
By Lemma [TJ we get 

\anK COs{y/\^t)Lpn{xQ)\ = '^\an\nfn{xo)\ 

= \^lnfn{xQ)\ 

< C£{AfY'^. 

According to the Weierstrass M-test, ^ Q!„A„ cos(v%j"t)(/3„(a;o) converges uniformly 
for any t. Similarly for the term |/3„VXi sin(v%i'0</'n(a;o)|- This implies ut and uu 
exist in the classical sense. □ 

Remark. In [H], Hu used the eigenvalue estimate A„ = 0(n") to estimate the term 
'^\oin'^nCOs{^/\it)ifn{xQ)\. That argument requires slightly stronger regularity 
condition. If our argument is used to replace all eigenvalue estimates in that paper, 
we could arrive the following result: 

Let / be a real- valued function on R satisfying F{r) < C(l + |rp) where F(r) = 
/q f{s)ds. If gi G domA with Agi G dom£, and 52 G domA, then the nonlinear 
wave equation with Dirichlet boundary condition 

utt{x, t) = Au{x, t) + f{u) {x G K\Vo, t G R) 

u{x,Q) = gi{x) {xeK) 

ut{x,0) = 52(2^) {x e K) 

u{x,t) =0 (x G T^o, i G R) 

admits a weak solution, where the second derivative of u is the Frechet derivative 
of Ut in L^. 
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4. Finite Difference Method 
The wave equation on Tm is defined by 

Umix,t + 1) = 2Um{x,t) - Um{x,t - 1) + h^fi:^^^^HmUm{x,t) 

where hm is the time span. In this section, we find the diflFerence between solutions 

of the wave equation on K and r„i. 

First of all, wc prove that the wave equation on the approximate graph is stable. 

Lemma 3. Let V be a finite dimension inner product space. Let H be a positive 
self-adjoint operator on V with eigenvalues < 3. Let £h{u) = {u,Hu), h he a 
function on F X N and g he a function on V. Let u he the solution of the wave 

equation 

u{t + 1) - 2u{t) + u{t - 1) = -Hu{t) + h{t) {t > 1) 
w(0) = 
u{l) = g 

Then we have SnHt))^/^ < 2(||<7|| +ELi M^W- 

Proof. Let {«„} be the orthonormal eigenvectors of H with corresponding eigen- 
values A„. 

For the case h = Q, let g{x) = ^ UnVn- Then the solution is 
u{x,t) = y^a„ sm{dnt)vn{x) 
where On = cos~^(l — ^) and q„ = a„/ sin(^^„). So, the energy at time t is 
£H{u{-,t)) = ^Q2^„sin2(^„t) 

- Z.sin(0„)2^" 

,2 



^ 

1-^ 



4 

2 _ /Ill„ll2 



By the assumption A„ < 3, so we have £H{u{-,t)) < 4^0^ = 4||£f| 

For the general case, let Wg{x,t) be the solution of this homogeneous equation 
at time t with initial velocity g. The result follows from the formula for the general 
solution: 



t 



u{x,t) = Wg{x,t) + y2'^h(;k){x,t - k). 



k=l 

□ 



Next, we estimate the difference between a finite energy function and its step 
function approximation. Recall that ip™{x) is the harmonic extension of the func- 



tion 5xp on V„, 



Lemma 4. For f € dom£ , we have 

E / l/(^)-/(2/)l''^y<«axCax£(/), 

where /Umax = max ft, r^ax = maxr, and K^ ^ = suppV'™. 



INFINITE PROPAGATION SPEED FOR WAVE SOLUTIONS ON SOME P.C.F. FRACTALS 6 

Proof. Using \f{y) - f{x)\^ < C£{f), we have 

\fiy)-fix)\^dy<C £{f). 



'K 

Applying the contraction mappings Fj^ on both sides, we get 



\foF-Hy)-foF-\x)\'dy = / \.f{y)-fi^rdy 

= Cfi^r^SifoF-^). 

Thus, for any finite energy function / with support in F^K , we have 

l/(y)-/W|'dy<CMSaxCaxf(/) 

where m — \uj\. For x Cz Fi^Vq d Vm, K^.x is contained in a m — f celL Thus, 
1/(2/) - f{^)\'dy < C f^"nl-^'rZ^^'£{f\K^J 



K 



Summing the inequahty over Vm, we have 



Since Km,x covers K at most N times, J^xev '^(/l-ft^x) — ^ ^(/)- '-' 

We define {u,v)m — '^xev ''^{^)''^i^)l^m,x- Under this inner product, the oper- 
ator h'^fJ'^x^rn is self-adjoint. 

Lemma 5. For any f e dom£, we have 

iii/iu-ii/iIlJ < c./ji^jz;j{r). 

Proof. By direct calculation, we have 

|ll/ll™-ll/llij 

^ J2 [ \fi^)-fiy)\\fi^) + fiy)\dy 
^ in l/(^)-/(y)Pd2/] (e/^ \f{^) + f{y)?dy 



1/2 / s 1/2 



Then the result follows from Lemma |4] and 

E / \f(.^) + fiy)\'dy < 2Ni\\f\\'y^ 

< 2iV(||/||v„ + ||/|UJ 



□ 
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Theorem 6. Assume f G domA with A/ £ dom£ and g G domA. Assume both f 
and g vanish on the boundary B. Assume B C Vm and eigenvalues of —h^fi^^^Hm 
are < 3. Let u,„ be the solution of the wave equation on r„j : 



{Urn{x,t+l) = 2Umix,t) - Umix,t - 1) + h'^^;;^]^H,nUmix,t) {x ^ Vrn\B) 

Um{x,t) =0 {X e B) 

Um{x,Q) = f{x) {xeV,n) 

Umix, 1) = f{x) + hg{x) + !YlJ-7n,xHmUm{x, 0) {x € Vm) 

Then, we have 

\umix, t) - u{x, ht)\ < Ct{h^ + Vm"?-™^) {xeVn,., t G N) 
where u is the solution of the wave equation on K . 

Proof. Assume g = for simplicity. Let u = q;„ cos(-\/A^t)v3„. By Theorem [5J 
the classical solution u exists and £{Au) < oo. The discrete wave equation on Vm 
comes from discretization of uu and A as follows: 



u{x, h{t + 1)) - 2u{x, ht) + u{x, h[t - 1)) 
« h'^uu{x,ht) 
= h^Au{x,ht) 

~ h^lJ-m^.xHmU{x, ht) 

So we want to estimate the error that appears in those two discretizations. 
For the first error, let 

erri {x, t) = u{x, h{t + 1)) — 2u[x, ht) + u(x, h(t — 1)) — h^Utt{x, ht). 

We have 

erri(a:,t) = 2^ (cos{\/Kih) - 1 + ]^\nh^ \ anCos{\^ht)(pn{x). 
Using I co%{\fXnh) — 1 + \\nh^\ < ^A^/i'*, we get 

llerrilla = ^ 4 cos(v/A^/i) - 1 + -An/i^ I cos( VA^/imn)^ 



n=0 

N 



f44 

n=0 n=N+l 



for any N. By [12, Thm 4.1.5], A„ = for some a > 0. So, we can choose 

A^v = ©(t^)- Thus, llerrilll = 0{h^). Similarly, we have £{erri) — 0{h'^). Using 
Lemma [SI we have 

||erri||„ = 0{h^ + yJJlF^h^). 
For the second error appears in A w ^.^^Hm, let 

err2(a;, ht) — Au{x, ht) — ii^^H„iu{x, ht). 
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Using HmU — / Au ^l™^d/ip^ A. 2. 5], we obtain 

I|err2||™ = ^ fi']^ {Au{x) ~ Au{y)) i^i"''> {y)dy fim,x 

< f \Auiy) - Auix)\' i^i"^Hy)dy 

< V / \Au{y)^ Au{x)\^ dy. 
Using Lemma m we have 



llerralU = 0(77^^^). 

Let e{x,t) = u{x,ht) — u„i{x,t). Then, e satisfies the graph wave equation: 

e{t + 1) - 2e{t) + e{t - 1) = /iVm!x^me(t) + erri(t) + /i^err2(t). 

Also, we have e(0) — and ||e(a;, l)||m — 0{h^ + y/ii"^r"^h'^) by similar estimates. 
Thus, Lemma [3] implies 

= 0{th'^ + t^/JjiPH^h). 
And the result follows from ||e||oo — 0(fm(e)^/^). □ 

Example. In Sierpinski Gasket with uniform measure, it is known that Ex- 
ample 3.7.3] 

2 degfa:) ^-^ 

Since A„i/ is a graph Laplacian, the eigenvalues of —Am are less than or equal to 2. 
Since the condition of Theorem [5] is satisfied for < 5^™, we take hm = 5^™/^. 
The difference equation becomes 

ujhmjn + 1)) " 2u{hmn) + u{hm{n - 1)) _ 3 „ 



Note that the constant y ^S™ is the scaled propagation speed. In [0, 1], the constant 
is 2™, which is the inverse of the grid size. Thus, the propagation speed in [0, 1] 
is same for all m but it increases as m increases in SG. And this gives a heuristic 
reason that the wave in SG doesn't have finite speed, which was first observed in 

El- 

5. Infinite Wave Propagation Speed And Heat Kernel Lower Bound 



In this section, we use heat kernel estimate and a relation between wave and 
heat equations to obtain some off-diagonal behaviors for the wave equation. Since 
we need Neumann heat kernel estimate, we assume i? = in this and next section. 
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Lemma 7. Assume f G domA with Af G dom£ and g G domA. Let u he the 
solution of the wave equation. Let v{x,t) — -^^exp{—j^)u{x, s)ds. Then v is 
the solution of heat equation: 

J vt{x,t) = Av{x,t) {xeK,t£R) 
\ «(x,0) = fix) (xeK) 

Proof. By theorem [51 the classical solution u exists. Since 



OO 

u = ai+ Pit+ 2_^an cos(V A„t)<y9„ + /3„ — — (y5„, 

n=2 n=2 V'^n 

the energy £{u{t),u{t)) + \\ut{t)\\2 < A + Bt^ for some A and B. Thus, ||u(t)||oo < 
A' + B't. Since N{t,s) — ^ rapidly as s — >■ oo, v is well-defined and the result 
follows by direct verification. □ 

In [15j, Adam Sikora proved that for a large class of self-adjoint operator with 
Gaussian off-diagonal estimate, the heat kernel estimates are related to the propa- 
gation speed of the wave equation. For homogeneous hierarchical fractals, we have 
sub-Gaussian estimatefT] and this makes the propagation speed infinite. 

Theorem 8. Suppose the heat kernel satisfies the sub-Gaussian lower bound: 

p{x,y,t) > C exp{~^) {x,y £ K, l> t> 0) 

where j3 < 1. Assume f G domA, Af G domS , f > 0, f =/= and g = 0. Let u{x,t) 
be the solution of the wave equation. Then, for all x £ K and S < I, there is t < S 
such that u(x,t) > 0. 

Proof. Let v{x^t) = ^^^o ^-''^P(~ft )"(^' '*)^'^ defined in Lemma[71 Since v is 
the solution of the heat equation with initial value /, for < < 1, we have 

v{x,t) = j p{x,y,t)f{y)dy 

> C j exp{-^)f{y)dy 

= C||/||iexp(-i). 

Take x be any point in K. Suppose, on the contrary, u{x,t) < for t < 1. Since 
g — 0, we have supj ||u(i)||oo < A for some A > and 



OO 

f I S^ 

v(x,t) < 2 exp( )u(x,s)ds 

s 

OO 

r 1 

< 2A / ^=exp( )ds 



s 



It leads to a contradiction that C||/||i exp(--^) < 2Aexp(-|^) + 0(t^)) for 
t < S because P < 1. □ 
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Figure 5.1. Functions in Example [9] 



However, the wave oscillates in space instead of spreading, as will be illustrated 
by the following example. And this says we cannot expect u{x,t) to be positive 
within short times even if / > 1 and <? = 0. 

Example 9. Consider the Laplacian with Neumann boundary condition on SG. 
Using spectral decimation|14[ l6]. we can have the estimate 

-1.5 < ^4 < 2, 

where tp^ is shown in Fig [51 Now, we define / by combining copies of 1^94 as shown. 
On each level, the solution of the wave equation is of the form coa(\/\t)ip. So, 
the wave oscillates faster on the upper level. Let f — Af + 7 and u be the wave 
equation with initial position /. The classical solution exists even though — 00. 
Although / > 1, {2 is not positive even in a short time interval because ip^ — 2 at 
some point. 



6. Wave Kernel And Heat Kernel Upper Bound 

Since the wave solution has infinite propagation speed for some fractals, we 
would like to get off-diagonal estimates of the solution of the wave equation for 
those fractals. 

We define Ptu be the heat solution with initial data u after time t where i > 0, 
that is, 

PtC^anifin) = y^Q„e~'^"V„ 

where u = ^ cin^n- Also, define WfU to be the solution of the wave equation with 
initial data u and initial velocity after time t where < > 0, that is. 



In this section, we assume the heat equation satisfies the following kernel upper 
bound: 



(6.1) 



for some a and some /3 > 2 which is true for many fractals [HI [7] . 
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Lemma 10. Assume the heat kernel satisfies the upper bound (|6.ip . For f e 

L^{K), we have 

\Pt+i/zf{x)\ < ^exp {^-Cr^z^Rez 
for t > where r — d{x, suppf). 

Proof. By scaling /, we may assume ||/||i = 1. Let u{z) — Pi^i/^f(x) which is 
analytic on {Rez > 0}. Note that 

p{x,y,z) = ^er^"^Lpn{x)(pn{y) 

< 



(Rez)" ' 

Therefore, we have \u{z)\ < On the other hand, the kernel upper bound tells 
us that 

\u{z)\ < ^exp (-Cr/^z/^ 

Because of symmetry, we only prove the statement for the first quadrant. Consider 
the strip fl ^ {x + iy : < y < ^} and let 

v{z) = u(exp(z)). 

By assumption, \v{z)\ < for z e and \v{x)\ < exp{—C exp{ax)) for x G R. 
Let 

^ exp(-Cr/^ exp(-g3Yz)) 

Now / is analytic on the strip 51, continuous on f2 and bounded by 1 on boundary 
of 57. Also, it satisfies a decay estimate 

|/(z)| <exp(Cr^exp(^|z|)) (z e f!). 

By the Phragmen-Lindelof theorem, we have |/| < 1 on 51. Thus, for z G {Rcz > 0}, 
we have 



C g 1 

Hz)\ < —\exp{-Crfi-^zfi-\ 



< 



^ ( 



■Crl^ z'^ sin(— — |argz| 



< — exp (^-Cr -9-1 2-3-1 Rez 



□ 



Recall that the heat equation is the averaged wave equation and we can use this 
to recover the lower frequency of the solution of wave equation. Therefore, mollified 
solutions is exponentially small outside the support of initial function when time is 
small. 
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Lemma 11. For -h > J > and~l < a < 0, we have 



mm 

Re(z)=i 



z-^ +t 



z + t 



where Ca is a constant depends on a only. 

Proof. We may assume arg( ^_j;^^ ) > because of symmetry. For < arg( ^_}^^ ) < 



■|, we have 



z-^ +t 



z-i +t 



> 



1 

71 



For arg( ^_} ) > let z = 7 + ix. We have 



1 



7^ + 



z-^+t (t(72 + a;2) +7)2 + 
Since arg( ^_j;_^^ ) > j, we have x > t{-y'^ + x"^) + 7. Hence, 



7^ + .T^ ^1 



(t(72 + + 7)2 + -1-2 



Therefore 
1 



2-1 +t 



z + t 



> 



> 



Q + l 



((^(7^ + x^) + 7)2 + x2)"/2(t(72 + a;2) _^ 
a;"(te2 +7) 



where the last line comes from minimizing x over a; > 0. Combining the two cases, 
we get 

a + l 



mm 

Re(z)=7 



1 



Z"^ +t 



1 

Re( — ^— — ) > CQmin( 

z + t 



1 



-f-^+t 
> Cat'"/2min(t-\7)i+"/2 



^1+q/2^-q/2-j 



□ 



Theorem 12. ylssitme f/ie /leai kernel satisfies the upper bound ()6.ip . Let 0o-(t) = 
-7=^ exp(— 1^). For a = O ( ( ■^t^- j ^ j and f ^ , we have 



\i^,^Wf){x,t)\ < 



C 



a+l + 



2 j-J-^u5/2 



r.2 \ /3-2 



exp(-C - a)||/||i 



where r — d{x, suppf ). Furthermore, for a = O 
we have 



y-pw J 



\Wf{x,t)\< 



C 



^"+1+^^^0-5/2 



exp(— C 
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^2 \ 3-2 



and P-af e , 
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where r = d{x, suppf). 

Proof. The relation between heat and wave equation can be written as 



exp(-sA„) = / coa{ty/Xn)e '^'dt. 



Changing some variables, we get 



The inverse Laplace transform implies 

^cos(vA„<) = TT- W-exp -— + si|ds. 

for any 7 > 0. The mollified cosine is 



(j) * cos{y/X^-)j (t) ~ e 2^"cos(V%iO 



t r+'°° /tt / , , 1 



7 — 200 



^exp( -A„(^ + |) + s<2 ) ds. 



Hence, the mollified solution can be computed as 



2ni 



7 — 200 



(4s)-i + - 



■Jds. 



Rewrite the equation by letting u{z) — P(^4z)-^+f f ^iid v(t) = {(p^ * M^/)(t), we 
have 



vit) 



t 



1 



7+ioo 1 
-ioo S ^ + cr 

t? J- —L r-r 



7 — 200 
7 + 200 





1 








+ cr 



2-/3 



Re- 



Using Lemma [TUl we have 

w( — \ ) < — 

Using Lemma [Til for i > 47 > , we have 

1 C ( /3 P g-2 

I )| < — exp -Cr'3-i72/3-2cr2/3-2 

s -'- + cr <" V 

for Res = 7. By the Cauchy integral formula, we have 

I" ( _i , )l < -^-r;; ^exp(-Cr/^75F=2CT2F^ 



S-l +cr'(i +ct)2 s2 



(is. 



< 



exp - 

i"7 V 
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for Res — 7. Substituting in ()6.2p . we get 



bWI < [ h + is\'^^^ds + - [ |7 + is|"^/^ds)exp(-Cr/^7^5^cr^) 

^" Jo 7 Jo 



Cq^^ 1 1 f f |9-2 



The first result follows from putting 7 = 2/^-2 • The second result 

follows from the identity 

{4> * cos(VA^-)) {t) = e"^^" cos(VA^O- 



□ 
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